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Summary 

An asymptotic analysis is done for flow through and 
around screens with small resistance coefficients. Both 
steady and oscillatory flows are considered, but only the 
case of a screen normal to the flow is treated. At second 
order in the asymptotic expansion the steady flow normal 
to the screen is nonuniform along the screen, due to 
components induced by the wake and by tangential drag. 
Therefore, the third-order pressure drop is nonuniform 
and the wake contains distributed vorticity, in addition to 
the vortex sheet along its boundary. The unsteady drag 
coefficient is found as a function of frequency. 

Introduction 

A two-dimensional screen of finite height, placed into 
an otherwise uniform flow, acts as a resistant surface so 
fluid can pass through the screen only if driven by a 
pressure drop. Fluid which passes around the screen will 
not suffer a resistive pressure drop so that, by Bernoulli’s 
equation and continuity of pressure at the edge of the 
screen, it must have a higher velocity than the fluid in the 
wake. Hence, ideally the wake behind a screen is bounded 
by a vortex sheet (fig. 1). In general, vorticity will also be 
distributed throughout the wake. 

The strength of the wake vortex sheet depends on the 
porosity of the screen, as modelled by screen resistance 
coefficients. In this report an analysis is done for screens 
with asymptotically small resistance. Only a uniform 
screen normal to the oncoming flow will be considered. It 
is found that the distributed wake vorticity is small, of 
third order, in the resistance coefficient. 

An analysis of flow around a very porous screen was 
made by Taylor (ref. 1). He represented the screen by a 
sheet of sources whose strengths were chosen to give the 
appropriate pressure drop across the screen. Because he 
ignored flow induced by wake vorticity and effects of the 
sources within the sheet upon each other, Taylor’s 
velocity field is formally accurate only to O(k), where k 
is the small resistance coefficient of the screen. (By 
formally accurate we mean his approximation agrees with 
an exact asymptotic expansion to O(k). However, 
Taylor did not expand his solution in powers of k. The 
form in which he expressed his drag coefficient agrees 
with experiment up to £ = 4.) 


Taylor’s distributed source analysis was extended by 
Koo and James (ref. 2). They introduced a wake in an ad 
hoc fashion and ignored the flow in front of the screen 
induced by wake vorticity. Also, the inviscid condition 
that vorticity be a function of the stream function 
generally is not met by their model. Their model also 
produces spurious jumps of the tangential velocity across 
the screen and of the pressure across the wake boundary. 
All in all, Koo and James’ model does not improve on the 
formal accuracy of Taylor’s solution; however, it 
provides a good fit to data on drag coefficients (Graham, 
ref. 3). The present report is concerned more with the 
formal solution of governing equations than with robust 
modelling. 

The Approximate Analysis section of this report 
describes an approach which agrees with the Asymptotic 
Analysis section (p. 4) to O(k) in the velocity field. The 
approximate solution could be obtained by summing an 
infinite subset of terms in the formal asymptotic 
expansion; in the Approximate Analysis section the 
solution is derived by omitting the tangential drag on the 
screen and by ignoring the dynamical effect of wake 
vorticity. Thus, this approximation is analogous to 
Taylor’s but, because it includes different higher order 
terms, it produces a different solution. 

In the Asymptotic Analysis section a formal 
asymptotic expansion is carried out to 0(k 2 ). At this 
order effects of tangential drag and of wake vorticity 
have entered the solution, and the distributed wake 
vorticity is found. 

In the Unsteady Flow section (p. 7) the drag coefficient 
of a screen in oscillatory uniform flow is calculated. 

Governing Equations 

We wish to solve the inviscid flow equations, subject to 
certain jump conditions which model the effect of a 
screen on the flow. Nondimensional variables will be 
used in which the fluid density, screen half-height, and 
upstream velocity are unity. The resistance coefficient k 
determines the pressure drop across the screen [P] 
through 

[P]= l -kU 2 (x = Q,\y\<V (1) 


l 




Figure 1. — Defining sketch of screen between two walls. 


(ref. 1) where U(x=0) is the velocity normal to the screen 
and the screen lies along x = 0, |y|<l (fig. 1). The 
tangential resistance coefficient B determines the jump in 
tangential velocity [K]: 

[V}=BV + (2a) 

where V + is the tangential velocity on the upstream side 
of the screen. 

Mass conservation across the screen requires that 


Approximate Analysis 

In the next section an asymptotic analysis for k-~t 
be done. Here a less formal approach is used, althou. 
still is valid onl^y when Ar«l; it agrees with 
asymptotic analysis to O(k). This approximate anal 
illustrates the general approach of this report. 

Substituting U = 1 + u, V = v, and H=P + 1/2 (w 2 + 
into equations (3) gives 


[U] = 0 (2b) 

These jump conditions must be satisfied by solutions to 
the momentum and continuity equations: 


ox dx dy 


dy dx dy 


( 3 ) 


dU dV n 
~x~ + ~x~ —0 
dx dy 


For flow confined by wind tunnel walls at \y\ = d, an 
additional boundary condition is K=0 on the walls. The 
upstream condition is (/— 1, K— 0 as x oo. 


dH 

dx 



dH 

dy 


dv 

= — +«w 
dx 


du dv n 

T + a- =° 

dx dy 


(‘ 


Here u = dv/dx- du/dy is the vorticity, which vanishes 
outside of the wake. If co is set identically to zero in 
equation (4), //satisfies Laplace’s equation. Let R be the 
harmonic conjugate to H so that H + iR is an analytic 
function in the complex plane cut along the screen. If B is 
set to zero, then [//] = [/*] = 1/2 kU 2 (0, \y \ < 1 ) and 
[/?] = [K] = 0 (see eq. (6)). Now, by the Plemelj formulas 
(ref. 4) and the boundary conditions to the problem, the 
analytic function H + iR is found to be 


2 


outside the wake and 


H+ iR = ~rr 1 U 2 (0,y'/i) cotanh {y' -z) 

loo/ L 4a 

I - 2 J dy ' 


-tanh — (y' + z)- 
4a 


(5) 


where z = x + iy. This satisfies H~0 as x oo and R = 0 

on [y\ = d. The latter is required because, from equation 
(4) with w = 0 and the Cauchy-Riemann relations, 

v = R (6) 

It turns out that (/(0,[y|< 1) is independent of y. If its 
value is denoted by U s , equation (5) becomes 


kU 2 

H + iR= 5 In 

4m 


Sinh 2d 


sinh — (z + 1) 
2d 


kUj 
4 d 


(7) 


where the branch cut of the logarithm is on the screen. 
Using equation (4) with w = 0 gives 

U=-H+l- k -^ f(y) (8) 

where f(y) is an arbitrary function of y, except for a 
possible jump across the screen. To satisfy equation (2b) 
and the upstream condition U=l, 



inside the wake. 

Expression (12a) agrees excellently with the data in 
figure 7 of Koo and James (ref. 2) where data extend to 
A: = 9. Expression (12b) agrees with the data in their figure 
8 up to A: = 3, after which it falls below the data. Thus, 
when the present approximation is extended beyond its 
range of validity it gives wake velocities which are lower 
than observed. (Indeed, expression (12b) becomes 
negative at large k.) Equation (11) with d= oo is plotted as 
C D versus k in figure 2, along with data transcribed from 
Graham (ref. 3). In cases where k was not measured, 
Graham used the formula k = (d~ 2 - 1) to relate k to the 
screen open area ratio Laws and Livesey (ref. 5) state 
that the formula A: = 0.52 (fi~ 2 - 1) is in better agreement 
with measurements; the latter formula was used in 
transcribing the data. Also shown in figure 2 are results 
from Koo and James’ model and a Pade approximant to 
the two-term asymptotic expansion (see eq. (53)). For k 
greater than about 4 the flow through a screen becomes 
unstable and a recirculating wake may begin to form 
(Laws and Livesey, ref. 5). There is no basis for the 
present analysis at these higher values of k. 

In the solution (8), the wake upper boundary is the line 
y = 1, x>0. In the following asymptotic analysis, to 
lowest order, the wake rises above pi a distance s(.x) 
given by 


C 6 x<o 

Ay)= Jo *> 0 , \y\>l ( 9 ) 

0 x>0, [y|<l 

Hence, U is discontinuous across the wake boundary 
(x>0,[y| = 1), which is a vortex sheebof strength 1/2 kU 2 . 

U s is found by evaluating equation (8) on the screen, 
where U= U s . Thus, 

, r _ -2 + 2V1 + k-(k/d) 

s ~ k-(k/d) ( ’ 

Equations (6) to (10) determine the flow field for given d 
and k. For instance, the drag coefficient for the screen is 

C D = 2[P] = kU 2 (11) 

and far downstream the velocity is 

, ku s m ^ 

1 + —J (12a) 


s{x)= [ V (,v',y= 1) dx' = 




kUj 

4-u 


Re 


In - 


sinh — x' 
2d 


sinh — (x' +2 i) 
2d 


dx' (13) 


Equations (6) and (7) were used to obtain V. When d~ °° 


s(x) 

kU 2 


-1 
87 r 


(m 


4 Tan “ 1 ( 

*2 + 4 2 


) 


(14) 


In figure 3 s(x)/kU 2 is plotted for several values of d. 
Also shown are experimental data measured in the NASA 
Lewis 20 by 30 inch ( Boy Scout ) wind tunnel. In the 
experiments d=2 and A: = 2.0 and 0.6. Although the 
experimental wake boundary was very sharp initially, it 
became diffuse further downstream: the data indicate the 
height above the screen of the centerline of the shear 
layer. The data were normalized by the theoretical kU 2 
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O Graham (ref. 3) 

Pade approximate (eq. (531) 

Koo and James (ref. 2) model 

Equation (11) 


00 O 


Resistance coefficient, k 


Figure 2. — Drag coefficient as function of resistance coefficient 


O k ■ 2, d * 2 
□ k-0.6, d = 2 



.4 .8 


Figure 3. — Wake boundary as computed from equation (13) for values of d indicated. 


found from equation (10). It is seen that when k< 2 
theory and experiment agree. The experiments were done 
in collaboration with Mr. D. McKinzie (NASA Lewis 
Research Center) in the initial phases of a study on shear 
flow instability. 

Asymptotic Analysis 

For this analysis we return to equations (1) to (3), and 
an expansion in powers of k is sought: for example, 


U=Uo + kU\ + k 2 U 2 ... 

Also, B is assumed of the form 

B = kB ] +k 2 B 2 + k 3 B 2 ... 

Only the unbounded case, d=°°, is treated 
At lowest order 

U 0 = 1, -P 0 = 0, K 0 = 0 


At order k , equations (1) to (3) become 


-dP i dU x 

dx dx 


-dP x dV x 
dy dx 


dU | dV, 
dx + dy 


= 0 


removes the singularity (i.e., the 5-function associated 
with dU x /dy) from the first of equations (17). Thus, is 
given by equation (14) with ku}= 1. 

Now, if we let 

(15) H 1 = P 2 +\{u\ + Vty 

then equations (17) can be written 


with [ L/j ] =0 = [V\], [P j] = l/2. Equations (15) are 
essentially the same as equations (4) with u = 0. Their 
solution is found by letting d-~ °° and setting kU 2 = 1 in 
equation (7): 


P, + ,'Q,.— In 



(16a) 


— dH 2 d(J 2 
dx dx 


-dH 1 _dV 2 1 df 
dy dx 2 1 dy 


(18) 


and with equations (6) and (8) 


These equations indicate that V 2 H 2 = 0 except for a jump 
across the vortex sheets of 


U, 


iV x = ~ In 
47T/ 



U6h; {// 2 )= ^',(A-,Ly!=l)sgn(y) 


Note that this velocity is logarithmically singular at the 
edges of the screen. Hence, the expansion breaks down in 
a small region around the edges, order e~ l /k . 

In equations (16) the screen wake lies in x>0, [yj<l, 
because f(y) is nonzero there (eq. (9)). This representation 
of the wake would result in nonuniformity of the asymp- 
totic expansion at the edge of the wake because df/dy = 
-sgn(y) 6(|y|-l). To avoid nonuniformity, we replace 
/O') by f{y - sgn(y) v(x)) with s=ks x + k 2 s 2 + kh 3 + 
so that the upper boundary of the wake lies ony = 1 +s(x). 
This does not affect the first order solution (16). 

At order k 2 , equations (1) to (3) become 


-dP 2 

dx 


dUj 

dx 


1 ds\ 

2 dx 


V 

dy 


+ U i 


dU x 

dx 


+ V, 


d_Uj 

dy 


- dP 2 
dy 


dV 2 

dx 




dVf 

dx 


+ V, 


dV x 

dy 


dU 2 
~dx + 


5^2 

dy 


= 0 



The angled brackets denote a jump across the sheets. 
U\(x,\y\ = 1) is to be understood as the average of the 
velocities just above and just below |y| = 1, so 

(// 2 )=-Q + ^ J P,(x,l)) sgn(y) (19) 

Again, it is convenient to introduce the harmonic 
conjugate to H 2 ,R 2 \ then from equations (18), V 2 = R 2 . 
The jump of R 2 across the vortex sheets is found by 
removing a nonuniformity in the equations at next order 
(analogously to what was done with eqs. (17)). Doing so 
gives 




ds\ ds 2 
dx + dx 


dV j „ dR x 

= K 2 + *l -gJT =R2 + s ' -ft 


on the upper vortex sheet (y= 1, x>0). The jump of this 
equation is 


Taking 


(R 2 )-^(U x )=\ v l {x,\)=^Q l {x,\) 


( 20 ) 


s\ = j 0 Vi(x',y= 1) dx' 


where ds x /dx- V x (x,\) has been used. Also on the lower 
sheet 
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On the screen conditions (1) and (2) give 

[tf 2 ] = (/,((>, [y | <1)=-J (21) 

by using equations (16) and 

[/? 2 ]=fl 1 Q 1 (0,ly|<l) (22) 

Since [P 2 \ = [H 2 \, to this order of approximation the 
pressure drag is uniformly spread over the screen. 

A sectionally analytic function with jumps (eqs. (19) to 
(22)) is (Roos, ref. 4) 


It follows from equation (25) and the preceding 
method of analysis that at third order the downstream 
velocity inside the wake is 


Uy(oo,y) = - [P 3 ] - U l (<x>,y)(J 2 (°°,y) 

= ^ - U 2 (0,y) (26) 

for t< 1. Hence, the velocity is nonuniform within the 
wake or, in other words, vorticity is distributed 
throughout the wake. 

The total third-order drag coefficient of the screen is 
Q> 3 = I *_ ^ 3 ] dy = ^ + \[ fJiM dy 


H 7 + 


iR 2 =_L -WiZ'Y 

2 m |^8 \ /z + 1 / 4 \z + i / 


1 f- P x (x',\)-iQ x (x\\) J , 
— r ; : dx 

47 rt Jo x -z~ 1 


+ j_ r p { (x\\)+iQ x (x\\) dx , 

4tt/ Jo x'-z + i 


+ 


B\ 

2ir 


f> iQ\(Q,y ) 

J - 1 iy ' - z 


(23) 


To satisfy equation (2b) and the upstream condition U 2 , 

y 2 ~ 0: 


j 7 - \ H 2 (x\0,y) dy 
16 - 1 


(27) 


Substituting equation (16a) into equation (23) shows that 
the integrals required to evaluate equation (27) can be 
done in closed form. The result is 


C 




4-&(- 4(l " 2)2+4 f as) 


71= 1 


Thus, to third order. 


kl ski 

C D = k - — + — -0.08335,4:3 
2 4o 


(29) 


U 2 + iV 2 = 


- H 2 - iRj + -/(y) 


(24) The net tangential drag on half the screen is 


In equation (23) the first logarithm has branch cuts on the 
vortex sheets and the second is cut on the screen. 

It follows from equation (23) that at this order of 
approximation U is not constant on the screen; hence, the 
approximate analysis of the previous section breaks 
down. Also, since 


[ ’ t7(0 ,T)1 V) dy = [ ' 1 *2[ y 2 ] +kHJ x [ V 2 ] + k^Vy] dy 
Jo Jo 

to third order. Using equation (2a) gives 
[V 3 ]=B 2 V l+ +B 1 V 2 + 


1 1 Using equations (23) and (24) for V 2+ gives 

[^ 3 ]= 2 ^l(°.Lvl < l) + t72(0,b'l<l)= 32 +U 2 (25) 

and U 2 is a function of y on the screen, at third order the 
pressure drag is no longer uniformly distributed. This 
y-dependence of U 2 is produced by wake vorticity and by This is the tangential drag on a screen adjacent to a plane 
the tangential drag on the screen. wall. 
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Unsteady Flow 

The present approach can be applied to unsteady flow 
through a screen. Turbulent flow through a screen was 
analyzed by Graham (ref. 3). Here we consider an easier 
problem in which the unsteadiness is due to uniform 
pulsations of the flow. The governing equations (3) now 
have dil/dt and dV/dt added to their right sides. The 
jump conditions (1) and (2) are assumed to remain valid 
at the frequency of interest. 

The pulsating flow will be taken as 

U 0 = 1 + a cos ut (31) 


Substituting this and equation (31) into equation (32) and 
equating coefficients of ct n e imut give 

~ ^0/77 + \,q- 1 + 0m- !,<?'*' 0 m <70 ^ = 0m<7 + (36) 

for m> 0, where primes denote differentiation with 
respect to x. For m = 0, 

00 q = ~ 'Poq'l'- 2<Al ,q-\ (37) 

and <t> mq can be found by solving equation (36) 
recursively in q. 

When <7 = 0, because <j> m ,_i=0, equation (36) becomes 


with a< 1. Then, at 0(k) 


~dPi 

dx 


du { at/, 
= Uo ^ + ~dT 


-3P i 
dy 


= U 0 


dVi dV x 
~dx + "dT 


dU, dV x 
+ ~dV 


~ 0m - 1,0“ 0mO0 ’ - 0mO + i m u<t>m0 
Let 

4>\x\tt) + iQ(j> = - 0 ' 

(32) where 4> ' — 0 as x oo; that is, 


<t>(x;Q) = e - ** 



e iax 'ip{x') dx' + gO>) 

OO 




(38) 


(39) 


(40) 


with [P X ] = \/2UI and U\, Kj— 0 as -oo. As in 
equation (16a) 


P X = U 2 0 



-Tan -1 


\+y) 


(33) 


with (Jq given by equation (31). It is convenient to rewrite 
this as 


The function g(y ) is determined by condition (26) and is 
zero outside the wake. Then a solution to equation (38) 
can be found as 

0mo(*) = E A ?<Mx;ru>) m>() I 

r= 1 \ 

/ (41) 

000 = - 0M J 


p 1 = 


E '"’"(l) 1 ” 


n = 0 


(34) Clearly, a\ = 1. Substituting equations (39) and (41) into 
equation (38) and equating the coefficients of and 
to zero yield 


where 0 is the term in brackets in equation (33), with the 
^-argument suppressed. Comparing equation (34) with 
equation (33) gives 


0 -mn-'Pmm 000“!, 010 = 001=2, 020=1 


and all other 0/s equal zero. 

A solution to equation (32) will be sought in the form 


A? = 


r 

(m — r) 


A 


m- 1 
r 


r<m 


m m - 1 

E ^r~ _ 0/770+ E A 'r 1 

7=1 7=1 


(42) 


The right side of this last equation is 0 when m = 2. 
Hence, it follows that 


U i = 


gimut 


) W E(5)^M < 35 > 


/ot\\m 
\2 


<7 = 0 


£ A™ = 0 m>2 

7= 1 


(43) 
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Using equations (42) and (43) gives 


A? = 


(m —r)\ 


(44) 


with a i = 1 and 

m - 1 „ 

1 r a r 

r= l 

which completes the solution with <7 = 0. 

When <7=1, equation (36) becomes 

— 0 m + 1,0 4>m - 1,1 = 4 " ilflu ( 45 ) 

for m >0. 0 m + ] 0 is given by equation (41). A solution to 
equation (45), using equation (37), is found readily in 
terms of 0(x;ft) and 30(x;ft)/3f 2. The latter, which is 
abbreviated to 30, satisfies 


30 ' + z'O 30 = — i<t> 


(46) 


As x— 30/0 — 00, so the expansion (35) is 
nonuniform in x. However, it is easily renormalized: 
equations (41) and (47) can be combined as 

<7 = 0 L 1 J 

m m + 1 

+ ^ A?<Kra) + ^ J] + 0(a 4) (49) 

r = 2 r= 1 

with the x and y dependence suppressed. This form for 
the first term on the right side removes the 
nonuniformity. 

The preceeding solution for (J\ enables one to calculate 
the unsteady drag on the screen to 0(k 2 ). 

The pressure jump at second order is 

[Pj\ = U 0 U x = £ (|)'"' l (0 m +f0 w -i 

m = -oo 

+ l 4> "' + ')\ x= o e " VWI (50) 


Assuming 
m + 1 

<t>m 1= X< BT<t>(x\u>) + C,„u> 30(x;a>) (47) 

r= 1 

and equating coefficients of ct>(ra i), 30, and 0 in equation 
(45) to zero give 


C - 1 

m (w-1)! 


by having used equations (31) and (35) with 0„, defined 
by equation (49). This pressure is not constant on the 
screen, so the nonuniform drag occurs at lower order 
than in the steady case and is not due to wake effects. 

The second-order drag is determined by integrating 
equation (50) over the screen: 

° 2= £ (f < 51 > 

m = -oo 

with 


For r> 1 , 


Bf= -rA” 


m + 1 


r^m 


)(48a) 


m- 1 


B%=~B% + i- Y.B? r=m,B\=-2,B\ = 2' 

r= 1 


For r = 1 and m > 1 


gill _ (- in - 


+ c, 


n + B'l' 
m — 1 


‘+/lj " +1 


(48b) 


The B'{' can be determined recursively. Solutions for 
higher values of q are found similarly. 


D? = j ’ 0 m (O,>> ' ) + ^ 0 m _ j(0 ,T ' ) + ^ 0m + l (0,T ') dy' 
From equation (40) with equation (33) for \p(x,y). 


0(0)* J 0(O,^;fi) = 

cosh 2 ft 


cosh 2ft - 1 . _ sinh 2ft 

— — (7 + In 2/ft) — 

27T/ft v 1 4ft 


r sinh 2ft , . 

[chi 2ft - 7 - In 2/ft j + sht ( 52 > 


+ 


2 7T /ft 


where shi and chi are the sinh and cosh integrals and y is 
Euler’s constant (Abramowitz and Stegun, rel. 6). The 
term Dj can be evaluated using equations (49) and (52) in 
equation (51). In particular, the term linear in a is 



1.00 I— 



1 

2 


0(oo) + 0(a ) 


To second order, D + kD\=kl-D-i- Since [Pi] = l/2 Uq, 
= Uq. In the steady case, the Pade approximant 


D = 


kD { 

1-(A :D 2 /D x ) 


k 


\ + -k 


(53) 


fits experimental data quite well (fig. 2). Thus, it seems 
sensible to write 


Dl 

a 


[>- 


k<j>(<j)) - 
2 


■ + 0(a) 


(54) 


for the unsteady drag at frequency w. In figure 4 the real 
and imaginary parts of D'/o; are plotted against to with 
k=l, and in figure 5, |£Haj is plotted for k <4 with 
to = 0.1, 1, and 10. 



(a) Real parts of D l / a . 

(b) Imaginary part of D'/a. 

Figure 4.— Unsteady drag coefficient as function of frequency to with 
k= 1. 



Figure 5. — Unsteady drag coefficient as function of k with frequency to = 0. 1 , 1, and 10. 
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Concluding Remarks 

A method for calculating steady or unsteady flow 
through very porous screens has been described. Good 
correspondence is found between theory and experiment 
by carrying out the asymptotic expansions to second 
order and using Pade approximants. Muramoto and 
Durbin have applied the present method to flow through 
inclined screens (paper submitted to AIAA 18th Fluid 
Dynamics and Plasma Dynamics and Lasers Conference, 
July 1985). In this case distributed wake vorticity appears 
at second order in the asymptotic expansion; the 
theoretical wake velocity profile is found to be in 
agreement with experimental data. 
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